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Abstract 

Qh 

An active stream of literature has followed up the idea of skew-elliptical densities 
initiated by Azzalini & Capitanio (1999). Their original formulation was based on a gen- 
eral lemma which is however of broader applicability than usually perceived. This note 
examines new directions of its use, and illustrates them with the construction of some 
probability distributions falling outside the family of the so-called skew-symmetric dens- 
ities. 
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1 Background and aims 

A currently active stream of literature deals with d-dimensional continuous probability distri- 
O butions such that their density function can be written in the form 

> 

f(z) = 2f (z)G{w(z)}, zeU d , (1) 

where f , called the 'base function' in this paper, is a density function satisfying some form of 
symmetry condition, and G and w are functions whose requirements will be recalled shortly. 
The more prominent representative of this formulation is the skew-normal distribution whose 
density function at z is 

2(p d {z;D.)<S>{a T z) (2) 

where cp d denotes the <i-dimensional normal density N^(0,Q) and <5 is the scalar standard 
normal distribution function (Azzalini & Dalla Valle, 1996). In the general formulation D. is 
a covariance matrix, but for the present purposes we can restrict ourselves to the case where 
Q is a correlation matrix; a is a vector of parameters which regulate the skewness of the 
distribution. Clearly, setting f {z) = <pd(z;D.), G = O, w{z) = a T z in (1) produces (2). 
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In fact the chronological order of development of the two expressions above was the 
opposite. Near the end of a paper dedicated to the properties of distribution (2), Azzalini & 
Capitanio (1999) delineated a more general formulation, starting from the following result, 
whose proof is reproduced here because of it is extremely simple, yet instructive. 

Proposition 1 Denote by G the distribution function of a continuous random variable whose 
density function is symmetric about and by Y = (Yi,..., Y d ) T a continuous random variable 
with density function f . If the real-valued transform w{Y) has symmetric density about 0, then 
(1 ) is a d-dimensional density function. 

Proof. If X ~ G, independent of Y, the distribution of X- w (F) is symmetric about 0, implying 
that 



The first formulation descending from the above proposition assumed fo to be an ellipt- 
ically contoured density centred at and w{y) a linear function, leading to what was later 
called the family of skew-elliptical densities. Although the ensuing discussion mentioned 
that fo does not need to be elliptical, the actual development of Azzalini & Capitanio (1999) 
focused on this case. The idea of skew- elliptical distributions has been followed up and ex- 
panded by a several authors, including Branco & Dey (2001), Genton & Loperfido (2005), 
and a number of contributors to the book edited by Genton (2004). 

In these developments, it emerged that many results could hold replacing the assumption 
of elliptical distribution for f by the weaker assumption of central symmetry. From Serfling 
(2006), recall that a random variable Y is said to be centrally symmetric about if -Y has the 
same distribution of Y; if the density fo, say, of Y exists, then foiy) = /o(-y). A quite general 
formulation is expressed in the following result presented by Azzalini & Capitanio (2003). 

Proposition 2 Denote by fo(-) the density function of a d-dimensional continuous random vari- 
able which is centrally symmetric about 0, and by G a scalar distribution function such that 
G(-x) = 1 - G{x)for all real x. If w{z) is a function from U d to U such that w{-z) = - w(z)for all 
z e U d , then (1) is a d-dimensional density function. 

An essentially equivalent result has been obtained independently by Wang et al. (2004), 
with the component G{w(y)} in (1) replaced by a function n{y) which must satisfy the con- 
ditions < n{y) < 1 and n{y) + n{-y) = 1 for all y e U d . The term skew-symmetric density has 
subsequently been adopted to denote this set of densities, or equivalently (1), when the base 
function / is centrally symmetric, to emphasize the broader settings which includes those 
with elliptical density as as subset. 

An important result concerning these distributions is the existence of a stochastic repres- 
entation for a random variable Z of this type given by 



where X and Y are independent random variable with distribution function G and density 
function / , respectively, and the functions G, fo and w satisfy the conditions of Proposition 2. 
This representation is important not only for random number generation, but also for ob- 
taining further theoretical conclusions. An especially important implication is a distributional 



\ = P{X < w(Y)} = E Y {P{X < w(Y)\Y}\ = 



[ G{w{y)}f {y)dy. 



(3) 




(4) 
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invariance property stating that, for any function r(-) from U d to M q such that t{y) = t{-y) for 
all y, then t{Z) and t{Y) have the same distribution, written as 

t{Y) = t{Z). (5) 

Much work has been dedicated in recent years to distributions constructed via Proposi- 
tion 2, or its counterpart based on the function n{x). In addition to those already quoted, the 
review paper of Azzalini (2005) provides many other references, and since then the list has 
increased substantially. 

However, Proposition 2 is a special case of Proposition 1. Its role can be viewed in identi- 
fying some very simple conditions which ensure the fulfilment of those of Proposition 1. The 
aim of the present note is to examine situations which fall under the setting of Proposition 1 
but not of Proposition 2. More specifically, we shall consider cases where /o is not symmetric 
about 0, or w (•) is not an odd function. As a side effect, this exploration leads to a deeper un- 
derstanding of what has been developed so far in connection to Proposition 2, and it produces 
a more general formulation of representation (4) and of the invariance property (5). 

Before tackling the specific target of this article, we note that the statement of Proposi- 
tion 1 is still valid under somewhat weaker assumpions, as follows. We can relax the assump- 
tion about absolute continuity of all distributions involved, and allow G or the distribution 
of w{Y) to be of discrete or of mixed type, provided the condition P{X- w{Y) < 0} = \ in (3) 
still holds. A sufficient condition to meet this requirement is that at least one of the random 
variables X and w{Y) is continuous. 

In our development, we shall however work mostly with the original formulation of Pro- 
position 1, with only a few remarks related to the above weaker assumptions. 



2 Some general facts 

In the framework of Proposition 1, the very argument of its proof justifies the representation 
of a random variable Z with the stated distribution as 

Z=Y if X < w(Y) (6) 

where the condition X < w(Y) is satisfied with probability \. It is seen that the density of Z 
is proportional to the integrand function of the last term in (3). There is an obvious sample 
selection mechanism which underlies the transformation of the distribution f into /. 

A natural question is whether a representation similar to (4) can hold. The question is 
of theoretical interest, but also of practical relevance, since in random number generation 
the use of (4) in place of (6) avoids the rejection of half of the Y samples. Notice that the 
non-rejection of samples in (4) is achieved by exploiting the symmetry of the distribution / . 

It appears that some additional conditions must be required to provide a solution to the 
above problem. One such set of conditions is as follows. If there exists an invertible trans- 
formation R(-) such that, for all y e U d , 

My) = folR(y)], |deti?'(y)| = 1, w[R{y)] = -w{y), (7) 
where R'{y) denotes the Jacobian matrix of the partial derivatives, then 

z \y ifx*wm, 

[R~ l m otherwise, 
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has distribution (1). In fact the density function of Z at z is 

fiz) = Mz)G{w{z)} + f {R{z))\detR'(y)\ [l -G{w{R{z))}] 
= Mz)G{w{z)} + f {z)[l-G{-w(z)}] 
= 2f {z)G{w{z)} 

using (7) and G{-x) = 1 - G(x). Under the weaker assumptions indicated at the end of Sec- 
tion 1, the latter equality may not hold for all x. However this may affect only a number of 
x values at most countable, and the density function /(•) can be replaced by a regularized 
version without affecting the distribution. 

According to the first two conditions in (7), the density fo is required to behave according 
to a "generalized symmetry" with respect to R{-), that is fl _1 (F) must have the same density 
fo of Y. In addition w{-) must be an odd function in this generalized sense, as required by 
the third condition (7). In Proposition 2, in (4) and in (5), the transformation function is 
R{z) = -z = R- l {z). 

If in (6) we reverse the sign of the inequality, this generates the dual variable of Z, having 
density function 

2/ (y)[l - G{w{y)}] =2f {y)G{-w{y)} = 2f {y)G{{w{R(y))} 

where the last equality makes use of the third condition (7) . 

For the analogous of the invariance property (5), consider a transformation f(-) from U d 
to U q which is even in the adopted generalized sense, that is 

t{z) = t{R-\z)) (9) 

for all z e U d . From representation (8), it is then immediate that (5) holds. 

The function R{-) which satisfy (7) does not need to be unique. Typically, if R is one such 
function, R' 1 is another one. In some cases, there are more than two functions R which 
satisfy (7), as we shall see in Section 3.1. Vice versa, in other cases i?" 1 = R, for instance when 
R{z) = -z, and we effectively have a single transformation. 

In the one-dimensional case a general way to construct a random variable with symmetric 
distribution is via its integral transform. Specifically, if Y has distribution function Fo, then 
F (F)-i is uniformly distributed on {-\,\) and application of Proposition 1 provides the 
following conclusion. 

Proposition 3 Assume that fo is a probability density function with distribution function Fo on 
So <= R, G is a distribution function over U which assigns a probability distribution symmetric 
about and it is continuous except possibly at 0, and w\ is an odd function on {-\, \). Then 

2f (x)G{w l [F (x)-l]}, 2/o(x)G{ M /i[i-F (x)]}, (xeS cR), (10) 

are density functions on So. 

The same result could however have been equally obtained from Proposition 2, although 
via a slightly more involved argument, as follows. Start by taking the base density function to 
be the uniform density in (-|, |), w{z) = u>i(z), and then shift the distribution to the interval 
(0,1); we arrive at the density 2G{u>i(x- ^)\. Finally, transform this distribution via F^ 1 (0, 
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to produce the first expression in (10). The second expression is obtained similarly using 

w{z) = Wi(-Z). 

For d > 1, the probability integral transform is not so easily tractable as for d = 1. A 
manageable case occurs with d = 2 and independent marginal components. If Y\ and Y 2 
are independent variables with distribution functions F\ and F 2 , respectively then a simple 
computation gives 

pit) = P{Fi ( Yi) F 2 (F 2 ) < r} = t (1 - log t) , (0 < t < 1) . 

Recall that p[Fi{Yi) F 2 {Y 2 )] is uniformly distributed in (0,1). If /1 =F[, f 2 = F' 2 and wi is an odd 
function, then from Proposition 1 

2/i(yi)/ 2 (y 2 ) G{m [piF^FAy^] - I)} 

is a proper probability density on the support given by the Cartesian product of the support 
sets of Fi and F 2 . 

This mechanism is of general validity for d = 2, but it has the disadvantage that the base 
function fiiyi) f 2 {y 2 ) does not allow for dependence between the marginal components. In 
addition, in many cases, the distributions so produced are not very convenient to work with 
from a mathematical viewpoint. For these reasons, in the following section we examine other 
mechanisms, which are simpler to handle. Moreover they are more directly related to the 
formulation expressed by (7)-(9); hence they better serve the purpose of illustrating the 
underlying mechanism. 

3 Some specific constructions 

We shall now examine a few specific distributions to illustrate the formulation of the previ- 
ous section. Given this aim, a complete study of the properties of these distributions is not 
attempted here. To ease exposition, we shall concentrate on the case d = 2. 

3.1 Cases with w{y) non-odd 

To start with a simple example, consider the probability distribution on IR 2 whose density 
function at z = {zi,z 2 ) J is 

f{z)=2(p 2 {z;Cl)<$>{a{z\-zl)} (11) 

where Q is a 2 x 2 correlation matrix with off-diagonal term p, and a is a real parameter. This 
f{z) bears a superficial resemblance to the skew-normal density (2) for d = 2, but with the 
noticeable difference from the skew-normal distribution that (11) is not skew at all; in fact 
it is centrally symmetric. The effect of modification of the normal density when a = 2 and 
p = 2/3 is shown graphically in the left panel of Figure 1. This density is bimodal, but other 
choices of a produce a unimodal density. 

The fact that the appropriate normalization constant in (11) is 2 cannot however fol- 
low from Proposition 2 which requires an odd function w, while w(z) = a{z\~ z^) is even. 
Equivalently n{z) = 4>{a(z 2 - does not satisfy the condition n{z) + n{-z) = 1 required in 
the formulation of Wang et al. (2004). However, if Y = (Y l ,Y 2 ) T ~ N 2 (0,n), it is true that 
a(F 2 - F 2 2 ) has a symmetric distribution about 0, and so Proposition 1 apply to conclude that 
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Figure 1: Contour level plot of density function (11) when a = 2 and p = 2/3 (left panel) and 
similar density when the last factor is ${zy - z 2 - {z\ - zf )} (right panel) 

(11) integrates to 1. In this respect, it would be irrelevant to replace <5 in (11) by some other 
distribution G with G' even. 

Much more general forms than (11) can however be handled, taking into account the 
following statement, whose proof is trivial and omitted. 

Proposition 4 If (X,Y) is a random variable on U 2 with density function h(x,y) which is a 
symmetric function of the variables, that is h{x,y) = h{y,x), and u>2{x,y) is a real-valued function 
such that u>2{y,x) = -u>2{x,y), then u>2{X, Y) has symmetric distribution around 0. 

A simple example of function w 2 which fulfils the above requirement is 

w{z) = w 2 {z 1 ,z 2 ) = a- l {z 1 -Z2) + --- + a m {z™-z 2 m ) (12) 

where z = {zi,z 2 ) T e M 2 and m is some natural number. The right-side panel of Figure 1 shows 
the contour level plot of the density obtained when the argument of $ in (11) is replaced 
by (12) with m = 2, a\ = \, <X2 = -1, and p = 2/3 as before. Notice that in general (12) is 
neither odd nor even, if the coefficients a\, . . . , a m are unrestricted. It is even if all a/s of odd 
order are 0, and it is odd in the dual case with only coefficients of odd order, we are back to 
the setting of Proposition 2. Combining Proposition 1 with Proposition 4, we can state the 
following corollary. 

Proposition 5 If fo(z) is a density function centrally symmetric over U 2 and fo is a symmetric 
function of the components of z, w 2 is as in Proposition 4 and G is a continuous distribution 
function over U which assigns a probability distribution symmetric about and it is continuous 
except possibly at 0, then 

2f {z)G{w 2 {z l ,z 2 )}, z={ Zl ,Z2) T eU 2 , (13) 



is a density function over U 2 . 



Figure 2: Contour level plot of density function (15) when a = \Jnl2 (left panel) and density 
function given by the first expression in (17) when p = -2/3 and a = \Jnl2 (right panel) 

For simplicity, we now restrict ourselves to the case where f {z) in (13) is ip 2 {z;Q.) and 
u>2 is as in Proposition 4, but the essence of the conclusions remains true for other centrally 
symmetric densities /o, symmetric of the arguments. The stochastic representation (8) for a 
random variable Z = (Zi,Z 2 ) T with such density function holds by choosing a transformation 
which swaps the two variables, that is 



R{y) = R y, 




and consequently also (9) applies for suitable f's. Among the many implications of this fact, 
the Mahalanobis distance Z T Q _1 Z is x 2 - Another consequence is that, since t{y) = yiy 2 = 
yiyi = t[Roy), then E{Zi Z 2 } = E{Yi Y 2 } = p. When central symmetry of f{z) holds, for instance 
when (12) is even, this implies that E{Z]} = E{Z2l = 0, and then cov{Zi,Z2} = p. 

We now move to a different example, and consider the bivariate distribution studied by 
Arnold et al. (2002) whose density function at z = {zi,z 2 ) T £ U 2 is 

f{z) = 2(p 2 {z;I 2 )®{az 1 z 2 ), (15) 

which is centrally symmetric too, and it enjoys several interesting properties. If Z = (Zi,Z2) T 
is a random variable with this distribution, then each of Z\, Z 2 is marginally standard normal, 
and the distribution of each component conditional on the other one is univariate skew- 
normal with parameter a, that is of type (2) with d = 1. Other properties are presented in 
the quoted paper, including an expression for coi{Z\,Z 2 }. A graphical representation of this 
density is shown in the left-side panel of Figure 2 for the case a = VtF/2 ~ 1.253 which has 
been proved by Arnold et al. (2002) to be the highest possible value producing a unimodal 
density; with larger a there are two modes. 

From the viewpoint of present context, (15) has a function w{z) = az\z 2 which is even, 
and again of Proposition 2 does not apply. It is however true that w (Y) has a symmetric 



distribution around 0, if Y ~ N 2 (0, I 2 ), and Proposition 1 confirms that 2 is the appropriate 
normalization constant in (15). More importantly, this explains that the normalizing 2 factor 
required in (15) is not by accident, and it connects the distribution of Arnold et al. (2002) 
with an apparently unrelated type of construction. 

To fulfil conditions (7) we can choose R{y) to be any of Rj-y for j = 1, . . . ,4 where 

«1 = (J "q 1 )- R2 = -Rl, ^3 = (" 1 °), i?4 = -i?3. (16) 

The first two of these matrices correspond to n/2 and -nil rotation, respectively, hence R 2 = 
Rj = R^ 1 and Rf = I 2 . This provides a stochastic representation for Z, and so a mechanism for 
random number generation. Using (9) with t{z) = z\+z\, we can state that Z 2 +Z| ~ j 2 . Using 
t{z) - {ziz 2 ) 2 we obtain that E{(ZiZ 2 ) 2 } = 1, which in turns implies cor{zf,zf} = 0, recalling 
that Zi and Z 2 have standardized marginals. 

If in (15) we want to replace I 2 by a correlation matrix D. with off-diagonal element p, 
then = a Y\ Y 2 does no longer give rise to a symmetric distribution around if (Yi, y 2 ) T ~ 
iV 2 (0,Q); hence Proposition 1 does not apply. The symmetry condition is fulfilled instead by 
w{Y) = aY\{Y 2 - pYi) and by the dual function w[Y) = aY 2 {Y\ -pY 2 ), since we have again the 
product of independent normal variables with mean. It then follows that 

2<p 2 {z;Q.)<&{az l (z 2 -pz l )}, 2<p 2 {z;Q.)®{az 2 (z l -pz 2 )} (17) 

are legitimate density functions on U 2 . The right-side panel of Figure 2 displays the density 
corresponding to p = -2/3 and a = s/n/2 as before. 

To search for a function R{-) fulfilling (7), we start by imposing the first and the third 
condition indicated by (7). For each given point zq e IR 2 , it is required to find a point which 
lies on the ellipse representing the locus of all points having the same density of zo and it also 
lies on the locus of all points having w (z) equal to -w{zo), that is to solve the equations 

D{z) = Dq, w{z) = -w 

where 

D(z) = z T DT l z = r(z 2 -2pziz 2 + z 2 ), w{z) = Z\Z 2 ~ pz\ 

1 — p 

and D = D{z ), wo = w{z ). After simple algebra one arrives at a fourth degree equation for 
zi with only even-order terms, or equivalently to the quadratic equation 

(l-p 2 )u 2 -D (l-p 2 )u+ u>l = 

where u-z\. The roots of this equation lead to four values for z\ - +\fu and four cor- 
responding z 2 = (_pz\ - wq)Iz\. There is a clear similarity with the earlier case where p = 
and four transformations associated to matrices (16), but now the four transformations R{-) 
are non-linear. The construction is illustrated by Figure 3 for the given point z = (2, 1) T , 
w{z) = z\ {z 2 - pzi) and p equal either to 1/3 or to 2/3; the change of pattern in the two cases 
depends on the sign of w(zq) as p varies. 

The final step is to check whether the second condition in (7) holds. Some algebraic work 
nor reported here shows that |deti?'| ^ 1. The implication is that, while (17) are valid density 
functions, a stochastic representation of type (8) does not hold for them. However (6) still 
holds. 



8 



Figure 3: Illustration of the possible choice of R[zo) for the first of the densities (17) when 
zq = (2, 1) T , a = 1 and p = 1/3 (left panel) or p = 2/3 (right panel). The ellipse represents the 
locus of points with the same density of zq; the dashed lines denote the loci of points with 
w{z) equal to -w{zq) 

3.2 Cases where / is not a symmetric density 

There are cases where Y is a continuous random variable with density function fo which is 
not symmetric about but there still are functions w{Y) with symmetric density about 0, so 
that Proposition 1 can be employed. Some constructions of this form have been sketched in 
Section 2. 

If d = 2, a simple but quite general formulation of this type can be obtained when /o(-) and 
w(-) satisfy the requirements for h and u>2 in Proposition 4, respectively. In other words, we 
are asking that 

fo(y) = fo(Roy), w{R y) = -w{y) 

where Ro is as in (14). Under this setting, (8) holds with R(y) = Roy. 

A simple illustrative example of this formulation can be started by choosing fo to be the 
product of two Gamma densities, such that each marginal density is 

/ 1 (x) = ^-x w - 1 e- x 
J r(w) 

if x > 0, and otherwise, for some o> > 0. For fo (y) = f\ {y\) f\ (j/2) where y = {yi, yi) 1 , there is a 
natural line of reflection given by the identity line. Among the many feasible functions w, we 
choose the simple linear case w{y) = a(yi -yi) where a is an arbitrary parameter. The density 
under consideration is then 

f{zi,z 2 ) = — -^(zi z 2 T- 1 exp(-zi - z 2 ) G{ec{zi - z 2 )} (18) 

for Zi,Z2 > 0, and otherwise. Here G is as required in Proposition 1 with possibly a discon- 
tinuity at 0. If Y = [Yi, Y2) 1 has density fo, then f{z\,Z2) is the distribution of Y conditionally 
on the event {X < a{Y\ - Y2)}, where X ~ G is independent of Y. 
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Various implications follow from (5) for a random variable Z = (Zi,Z 2 ) T having density 
function (18). Some examples are 

E{Z L + Z 2 } = E{Y y + Y 2 ] = 2oj, E{Z L Z 2 } = E{Y\ Y 2 } = a) 2 , 

irrespectively of G, since both t{y) = y\ + y 2 and t{y) = y\ y 2 satisfy (9). 

To compute the marginal distributions of Z\ and Z 2 , assume for the moment that a > 0, 
and rewrite the event {X < a{Y\ - Y 2 )} as {T < Yi\ where T = a~ 1 X+ Y 2 . Denote by F T the 
distribution function of T, which is the integral of the density 

poo 

Mt)= My 2 )g[a{t-y 2 )]ady 2 (19) 
Jo 

where g = G'. It is simple to obtain the density functions of Zi and Z 2 which are 

2/i(zi)F r (zi), 2/ 1 (z 2 ){l-F T (z 2 )}, (zi > 0, z 2 >0), (20) 

respectively taking into account that P{T < YiJ = \. When a < 0, the density function of T is 
still computed from (19), with a replaced by its absolute value, since the distribution of a~ l X 
does not depend on the sign of a, and the above expressions of the densities of Z\ and Z 2 are 
exchanged. Another way of looking at this aspect is to take into account (8) and to notice 
that reversing the sign of a corresponds to swap the distribution of the two components. 

The explicit computation of (19) is feasible for a suitable form of g. A convenient option 
is to set g equal to the Laplace density 

g(x) = i<r |x| , del), 

leading to 

frit) = 



e at ifr<0, 



2(1 + ar 

(e~ at h + e at I 2 ) otherwise, 



where, for integer a), 

Jo 



a) 



if a = 1, 



i tti-1 f k 



rM e -^-a)y 1 I ifa#1 

\ (1 - aT to *! (1 " uT- k J 

f*QO 



= T{a>)e- ta+a) 



tm\ + ar- k ' 

It is possible to produce a general expression of the distribution function Fj. However, for 
simplicity, we now restrict ourselves to the case cd = 1, such that 

1 



F T W = { 



2(1 + a) 



e at iff<0, 



l-e" r (| + |) if f>0, a = l, 

a 2 1 

1 - e~' + e~ at if t > 0, ajtl, 

a 2 -I 2(a-l) 
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Figure 4: Density function of Z in the exponential-Laplace case when w{y) = 2(y\ - y 2 ) in the 
left panel, and marginal density of Z\ (continuous line) and Z2 (dashed line) in the right 
panel 



when a > 0. By inserting this expression of Frit) in (20), we obtain the marginal density 
functions of Z\ and Z 2 . Figure 4 displays the bivariate density function of Z and the marginal 
density of Z\ and Z2 when a = 2. 

The moments of the marginal densities can be computed by direct integration, which for 
a > and any r > gives 



T(r + 1) 



r + A\ 
1 



E{Z[} = < 
E{Z 2 r } = 2r(r + l)-E{Z[} 



a 
~2J 



1 



a + ay 



if a = 1, 

if a > 0, a # 1, 



while the moments for a < are obtained by swapping the subscripts of Z\ and Z2 in the 
above expressions and taking the absolute value of a. From marginal moments up to second 
order and from the fact obtained earlier that E{Zi Z2} = 1, we can compute the correlation. 
Numerical evaluation indicates that the correlation increases monotonically from when 
a = to l/\/5 ~ 0.447 when a = ±00. 

We close this section with some remarks on possible extensions. The above discussion has 
focused on the case with linear w(_y) = a{y\ -y-i) but many other options are possible, such 
as w(y) = a{yf - y|) or w(y) = sin[a(yi - j^)]. Another direction is to adopt a base density 
with dependent components, since our choice of using independent components was only 
for mathematical simplicity. The overview of bivariate exponential distributions presented in 
Section 47.2 of Kotz et al. (2000) includes several distributions which meet this requirement. 
Among them, two especially important proposals are the Gumbel and the Marshall-Olkin 
bivariate exponential distributions, whose joint survival functions are 



exp{-y! -y 2 - Ayi y 2 ] , exp{-yi - y 2 - Amax(j/i, y 2 )} 



respectively, for y lt y 2 > 0, when scale factors are not included; here A is a positive parameter, 
and A < 1 for the first case. 

The Marshall & Olkin distribution presents the somewhat peculiar feature of a positive 
probability mass assigned to the line yi = y 2 . This situation is not handled directly by Proposi- 
tion 1, not even its extension stated at the end of Section 1, because here Y is not a continuous 
random variable. The argument of Proposition 1 can be applied to fo, if this denotes the den- 
sity function on the non-singularity set. The singularity set would not be affected, and a 
separate mechanism can possibly be adopted to modify the probability distribution on this 
set. 

4 Final remarks 

We conclude with some remarks on cases with dimension d > 2, when the base function fo is 
centrally symmetric. To ensure that (1) is a proper density we need that w(Y) is symmetrically 
distributed around when Y has density fo. The case when w (■) is odd, which of course 
includes linear functions, is covered by Proposition 2, and symmetry of w(Y) holds in general. 

For a non-odd w which is a function of two components of Y only, we can still make use 
of Proposition 5 with a d-dimensional fo, by using the following argument. If / is centrally 
symmetric, the same is true for the distribution of any subset of its components. Hence a 
function w of the nominated components of the form required by Proposition 4 produces a 
symmetric w[Y). 

We are then left with the case where d > 2 and w (F) is a non-odd function of at least three 
components of Y. The characterization of the set of such functions if is a problem which 
does not appear amenable to a simple solution, and it deserves a separate study, given also 
its interest as a property of centrally symmetric distributions. 
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